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We show that the one-loop effective action at finite temperature for a scalar field with quartic 
interaction has the same renormalized expression as at zero temperature if written in terms of 
a certain classical field <j> c , and if we trade free propagators at zero temperature for their finite- 
temperature counterparts. The result follows if we write the partition function as an integral over 
field eigenstates (boundary fields) of the density matrix element in the functional Schrodinger field- 
representation, and perform a semiclassical expansion in two steps: first, we integrate around the 
saddle-point for fixed boundary fields, which is the classical field 4> c , a functional of the boundary 
fields; then, we perform a saddle-point integration over the boundary fields, whose correlations 
characterize the thermal properties of the system. This procedure provides a dimensionally-reduced 
effective theory for the thermal system. We calculate the two-point correlation as an example. 



I. INTRODUCTION 

Effective actions are generating functional for the one-particle irreducible vertices of a field theory. From the 
vertices, one obtains all the correlations of the theory, which may be used to compute physical quantities, after 
suitable renormalization. In finite-temperature field theories [1] , the effective action is the thermal equilibrium Gibbs 
free energy of the system of relativistic quanta of the theory. Effective actions are, therefore, natural quantities to 
compute at finite temperature. For instance, effective action calculations [2-7] are among the techniques used to 
improve the poor convergence of the perturbative series in the bosonic sector, plagued by infrared divergencies. 

In this paper we build an effective field theory for quantum statistical correlations using the boundary fields as 
the relevant degrees of freedom, which we denote by boundary effective theory (BET). To do so, we start from the 
partition function of our field theory in equilibrium at temperature T = 1//3 written as a functional integral of the 
diagonal element of the density matrix (the Boltzmann operator for the field theory Hamiltonian) in the Schrodinger 
field-representation. The integral is performed over the eigenvalues of the field operator in that representation, which 
are the stochastic variables of the problem, and which we will henceforth call boundary fields. 

The density matrix element itself may also be written as a functional integral, one that describes an Euclidean time 
r evolution from a given boundary field configuration at r = 0, to that same configuration at r = /3. The integral 
described in the previous paragraph is then performed over boundary fields. The quantum statistical correlations of 
the boundary fields should fully characterize thermal equilibrium. 

The use of the density matrix functional leads to the derivation of the effective quantum statistical description for an 
underlying field theory, given by the effective action constructed from its microscopic Hamiltonian. The correlations 
of the stochastic fields (boundary fields, which only depend on spatial coordinates) express the thermalization encoded 
in the density matrix. 

The density matrix functional was already used in a similar context in [8], where it led to the construction of 
dimensionally-reduced effective actions. Likewise, in [9], it was used up to one- loop order to investigate the thermody- 
namics of scalar fields. In the latter reference, the boundary fields (considered as fluctuations around a homogeneous 
background) played a very nontrivial role in the calculation of the partition function of scalar fields. In quantum 
mechanics, this formalism was developed in Refs. [10-12] and led to the functional density matrix formulation of 
quantum statistics in Ref. [13]. 
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Our goal, in the present article, is to use the density matrix functional in a scalar field theory with quartic interaction 
to demonstrate a simple relationship between effective actions at zero and finite temperature, computed up to one- 
loop order. The relation provides a natural bridge between correlations, as well as renormalization conditions, at zero 
and finite temperature. In practice, it yields a recipe to read off the renormalized finite temperature result from its 
renormalized zero-temperature counterpart. 

The effective action that we obtain at finite temperature is a functional of the expectation value of the boundary 
(stochastic) fields. It is the generator of the one-particle irreducible vertices for the stochastic fields, and therefore 
allows us to reconstruct all stochastic correlations which define physical quantities in thermal equilibrium. In practice, 
our relation establishes a connection among those and zero-temperature physical quantities, such as particle masses 
and couplings. 

This article is structured as follows: in Sec. II, we write both the density matrix and the partition function as 
functional integrals, define the generating functionals for quantum statistical correlations, and associate them to 
thermodynamic free energies; in Sec. Ill, we integrate over the dynamical fields of the underlying theory, keeping the 
boundary fields fixed; in Sec. IV, we integrate over the (stochastic) boundary fields to obtain the effective action; in 
Sec. V, we go through the renormalization procedure to compute one-loop renormalized correlations; finally, in Sec. 
VI, we present our conclusions. 



II. DENSITY MATRIX FOR SCALAR THEORIES 



In quantum statistical mechanics, the partition function for a system in contact with a thermal reservoir at tem- 
perature T ((3 = 1 /T) is expressed as a sum (integral) over a stochastic variable, whose probabilistic weight is given 
by the density matrix. For a system described by a (self-interacting) scalar field theory, the stochastic variable is the 
field eigenvalue in the functional Schrdinger field-representation, 4>o(x), which satisfies <p\<j)o(x)) = (f> (x)\(j) (x)). 

The partition function is then a functional integral over the field eigenvalue of the diagonal element of the density 
matrix functional pp, 

z(j3)= [[VM*)]pii[M*),M*)], (!) 



Pp[Mx),Mx)] = (M*)\eM-PH)\<j> (*)) • (2) 

The Hamiltonian in the preceding formula specifies the underlying dynamics. In the present case, it is the dynamics 
of a (self-interacting) scalar field <j)(r,x). 

As is well known [14], pp can also be expressed as a functional integral over dynamical fields 4>(t,x) defined for 
Euclidean time r, subject to the boundary conditions 0(0, x) — <j>(/3,x) = 4>o{x), 

P p[Mx),Mx)]= J [mr,x)}e- s ^*)] . (3) 

4>(o,x)=4>(0,x)=4> o {x) 

The boundary conditions relate the stochastic variable of the integral in Eq. (1) to the dynamical fields that are 
integrated over in Eq. (3). The Hamiltonian 

H[U, <f>] = in 2 + i(V0) 2 + ^ <f + U(4>) , (4) 

which involves the time-dependent field <j>(r,x) and the conjugate momentum II(t, a:), leads to the Euclidean action 
(we assume d spatial dimensions), 

s[4>] = / (d D x) E m , (5) 

Jo 



m = \{dr4>? + \ W) 2 + ^ 2 + U(<j>) . (6) 

In this paper, we use the shorthands x = (r, x) and j(d D x) E = dr J d d x, with D = d + 1. 
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The density matrix is a functional of <fio(x) only. We may write 

pp[Mx),M*)] = e - 5 * wl , (7) 

Sd being a certain temperature-dependent dimensionally-reduced action. The field <fio( x ), the argument of Sd, depends 
only on the d spatial coordinates; all the r dependence of the original (d + l)-theory has been eliminated through the 
</> integration. The remaining integral over (f>o(x), required to obtain the partition function, is unrestricted (except 
for the vacuum boundary conditions that must be imposed at spatial infinity). 

The fields 4>o (x) are the natural degrees of freedom of the reduced theory. Any thermal observable can be constructed 
by integrating the appropriate functional of (j)o(x) over the fields (f>o(x) weighed with the corresponding diagonal 
element of the density matrix. The Euclidean time evolution can be viewed as an intermediate step which calculates 
the weights. 

Before proceeding, notice that the density matrix is not, in general, of the form e~^ Hd with Hd being independent 
of (3 as in ordinary quantum statistical mechanics; its /? dependence is far more complicated. This had already been 
pointed out in the analogous discussion of the transfer matrix carried out in Ref. [15]. The density matrix provides 
a direct but alternative way of deriving a dimensionally-reduced theory. 

In order to construct generating functionals for the dimensionally-reduced theory, we couple the field <j>o(x) to an 
external current j(x), obtaining a modified action, 



I[<h,j] = Sd{M - P J d d xj(x)M*) ■ (8) 

The associated partition function, 

Z\3\ = J {DM*)} (9) 

is the generating functional of quantum statistical correlation functions. 

The quantum statistical connected correlation functions are obtained as functional derivatives of the (Helmholtz) 
free energy functional F[j(x)} 7 defined as 

F[j] = -\ lim \ogZ\j] . (10) 
In particular, the expectation value of the field (f>o is given by 

<*,(*», = -g| • (id 

The label j in (<j)o)j is to stress that such an expectation value is a response of the system to an external current. In 
the sequel, we will drop the index j to simplify the notation. 

Another important quantity in the theory is the effective action T, 

T[{Mx))} - F\j(x)\ + [ d d xj(x) (Mx)) , (12) 



the Legendre transform of F[j\. The argument of the effective action is the expectation value of the field, an in- 
trinsic characteristic of the system. T is the generating functional of the one-particle irreducible quantum statistical 
correlations. We shall use it to impose renormalization conditions that connect with the physical (zero-temperature) 
parameters of the underlying field theory, and to implement a renormalization procedure. From a thermodynamical 
viewpoint, this effective action is the Gibbs free energy functional, essentially the pressure of the system as a func- 
tion of {4>o(x)). In magnetic systems, for instance, the effective action gives the dependence of the pressure on the 
magnetization. 



III. FLUCTUATIONS AT FIXED BOUNDARY 



The standard one-loop computation of the effective action for a position dependent background is performed in 
detail in Ref. [16] (Appendix 6 — 1). It makes use of the steepest-descent method, and can be thought as the first 
term of a semiclassical series. In the present case, one has a double integration to perform, which will require an 
adaptation of the standard techniques. 
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The first integral over <f>(r, x), in Eq. (3), will be dominated by configurations in the vicinity of classical solutions 
c satisfying 



□ E </> c + m 2 <f> c + U'(cj> c ) = 
(j) c (0,x) = <j> c (P,x) = <j) (x) 



(13a) 
(13b) 



where D E = — (d 2 + V 2 ) is the Euclidean D'Alcmbertian operator. Let us, for simplicity, assume that U is a single-well 
potential. Then, a unique solution 4> c (t, x) will exist for each boundary configuration <po{x). The classical solution is 
a functional c [0o] of the boundary field. We write 



(t,x) = 4> c {t,x) +t](t,x) , 
»7(0,x)=»7(/3,sb)=0, 



(14) 
(15) 



and expand the action around 4> c to quadratic order in r], 

5[0] = StfJ + S^S{<t> c , 77] + 6^S[cj> c , rj\ + <D{rf) 
The first-order variation S^S is given by 

6&S = J (d D x) E X - pMx) + m 2 <j> c {x) + U'{(j> c {x))\ V (x) + J d d x [ V (x) d T ^ c {x)f , (16) 

where we used the notation L4(t)]q = A(j3) — A(Q). The integrand in the first term of the right hand side of (16) 
vanishes identically because <j) c (x) obeys (13). The second term - a boundary term - vanishes due to the property 
(14) of the fluctuation 77. Therefore, <j> c is an cxtrcmum of S for ry-like variations (fixed boundary). 
The second-order variation 5^ S is 

5< 2 >S = \ J (d D x) E [[d, V (x))^ V (x))+ml V 2 (x) + U^M^)vH^) 

= \ J(d D x) E d,[r 1 (x)d^(x) 



+ 2 / ( d ° x ) E n B +m 2 + U"{4> c {x)) V (x) . 



(17) 



Using again (14), we obtain, up to quadratic order in 77, 

S[4>] « S[<f> c ] + \ J (d D x) E n{x) [d e + m 2 + U"(<f> c (x)j\ v(x) . 

It is convenient to introduce the Green's function 

[D E + m 2 + U" (&(*))] G[4> c ]{x; x') = 5^ (x - x') 
G[ ( f> c }(T,x;0,x , ) = G[ ( p c }(r,x;(3,x')=0 . 

In terms of G[0 C ], we obtain 

J [D4>{t, x)\ e- sw « e~ s ^ (det G[0 C ]) 1/2 . 

0(O,x)=0(/3,x)=0 o (a!) 

For single-well potentials it can be shown that S^S > 0, a necessary condition for </> c [0o] to be a minimum of S 



(18) 



(19a) 
(19b) 



(20) 



A. Building the classical solution 

In Appendix A, we obtain the following recursive relation for the classical solution </> c [</>o] 

f r -1/3 

4> c (t,x)= drx' <j> (x') d T 'G {T,x;T',x') 



- dr 1 J d d x' G (t, x; t', x') U' 



(21) 
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where Go(x,x') is a Green's function of the free operator 

p B +m§] G {t,x;t',x') = 5^\x-x') , (22a) 
G (r,x;0,a;') = G (t,x; j3,x') = . (22b) 

Diagrammatically, the classical solution 4> c {x) can be represented as the sum of all the tree diagrams terminated by 
the boundary field 4>o(x). It is also demonstrated in Appendix A that 

6 ^ T ' X) -[d T ,G[<j> c ]{T,x;T>,y)]l . (23) 



SMv) 



From the previous equation, if one knows the classical solution (j) c associated to some boundary configuration O , the 
approximate classical solution associated to the following fluctuation of the boundary condition, 

to(x)=^ (a:) + £(x) (24) 

is given by 

MMW = 0c(*) + / d d y£(y) [d^G^x-y , v )f Q + 0(f) , (25) 

as found in [9] . Another useful relation is 

^ J. — v- = [d T d T iG[(j> c \{T,x;T ,y)Y Q . (26) 

Finally, notice that 4> c [(j)o} is, in general, a non-periodic function of r. As long as the only condition on the 
configurations that enter in the calculation of Z[j] is to coincide at r — and /?, non-periodic configurations are 
allowed. 



IV. FLUCTUATING THE BOUNDARY FIELD 

From the first functional integration (see Eq. (20)) comes out an approximation for Sd defined in (7), 

Sdlfo] = S[MM + \ TrlogG-^cfo,]] ■ (27) 

Sd works as the classical action of the theory. Using again the recipe to obtain the one-loop quantum effective action 
(now for the expectation value of the field <po(x)), we arrive at 

0r[{<h)] = S d [(<h)] + \ Trlog<5( 2 )^[(0 o )] , (28) 

where Sd[(4>o)] is a shorthand for the second-order variation of Sd around (4>o)- However, to be consistent with 
the one-loop approximation, S^Sd[((t>o)] nas to be replaced by 5^S[(<j) }}. In fact, the expression for T as given in 
(28) contains only part of the higher-order terms. Therefore, we are left with the following one-loop effective action, 

pT[{<p )\ = S[<t> c ] + \ TrlogG" 1 ^] + \ Trlog 

where <j) c in this case is c [(0o)] • 

It is convenient to integrate the classical action by parts, obtaining 

S[4> c ] = J (d D x) E Q&D^c + \m 2 <fc + U(<f> c )j +lj ddx M*) [drMr, x)\o • (30) 

The presence of a boundary term is a direct consequence of the non-periodicity of 4> c [4> ] , pointed out in the previous 
section. Such boundary contribution introduces nonstandard terms in the calculations. 



5 2 S[<j> c 



(29) 
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At one-loop order, the expectation value of <f) is given by the saddle-point of the modified action I[4>o,j] (see Eq. 
(8)). Let us call <f> such an extremal configuration. Fluctuating the boundary field as 



o (x) = <j> (x)+£(x) , 
and denoting C [</><)] by <j> c , the first-order variation of I[(f>o,j] is given by 

5 W I[<t>o,j] = J (d D x) E pJ c (x)+m 2 Mx) + U'(Ux))] (fax) -$ c (x)) 

+ J d d xt(x) {[d T $ c {T,x)t - j{x)) . 



(31) 



(32) 



By construction, <j) c [4>o} satisfies the equation of motion (13), so that the first term on the r.h.s. of (32) vanishes. 
From (32), we obtain the saddle-point condition 1 for <p : 



(33) 



Notice that Eq. (33) is a condition on the time derivative of 4> c . Therefore, whenever j(x) ^ (even for constant 
currents) the classical configuration 4>c[4>o\ Wl U depend on r (even in the free case). 
For the second-order variation, we obtain 



where 



5 2 i[<f> ] = s 2 i (a) [M + s 2 i (b) [M , 

5 2 I {a) [M = \f (d D x) E (d D xy E (d D z) E [d Tl G^ c ]{z-y,x)]^{x) [n E +ml + U"Q c {z))} 

x [d T >G$ c ] (z; t', x')] q £(je') 



and 



5 2 I (b) [M = \ J d d xd d x'^x) [d T d T ,G[4> c ]{T,x;T' ,x')]l 

However, only the term 5 2 I^ [<f>o] contributes, because 

p E +m 2 + U"{4> c {z))] [d T ,G^ c ](z;r',x')] o =O. 
Using that the second-order variation of S[</> c [(fo]] and I[(f>o,j} coincide, we can write 

s 2 s\MM 



5cfro(x)5<fro(x') 



[drdr.Gi^X-y ,X% . 



We now have: 



pT[<j> ] = Sm + \ TrlogG- 1 ^] + \ Trlog ([d T d T ,G$ c f ) 
Using Eq. (19) written in terms of Go (sec Eq. (22)), we are led to 

[d T d T ,G$ c ]f = [d T d T , [g (1 + GoW'^y 1 }]^ . 
We may sum the series to obtain 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



[d T d T ,G[4>X = [drdr'G^ + [9 t Go]q U"{<j> c )G^G^ c ] [d T ,G f . 



(41) 



When j = 0, condition (33) finally implies that <fi c is a periodic function of r with period /3 (in fact, the trivial solution), but this 
property is not shared by classical solutions in the presence of nontrivial currents. 
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When we take the logarithm, we will have 

log [d T d T ,G$ c f = log So + log [1 + Go 1 (d T G U"@ c )Go 1 G\$ c ]d T ,G )] , (42) 

where we have denoted [<9 T <9 r 'Go]g as Go. Taking the trace and using that Tr (AB) n = Tr (BA) n to reorganize the 
series, we obtain 

Trlog [d T d T ,G\4> c ]\l = TrlogSo + Trlog [l + G^G\4> C ] {d T .G Q^d T G ) U"$ c )] . (43) 

We may write 

log (G" 1 ^]) = log (Go 1 ) + log (G-^JGo) . (44) 

Therefore, 

Trlog(G- 1 [^])+ Trlog (\d T d T .G$ c ]]l) = Trlog (G 1 )+ Trlo g g + 

+ Trlog {1 + (Go + d T G g 1 d T ,G )U"(4 c )} . (45) 

In Appendix B we show that the combination which appears multiplying U"(<j> c ) in (45) is the usual thermal propagator 
Af? which is given, in spatial Fourier space, by 

A F (k; r, t') = J- f(l + n{u k ))e-^- T '\ + n^e^^'A , (46) 
Zu k L J 

where Lu k = \J k 2 + TOq and n(uj k ) is the Bosc-Einstein distribution. In terms of A F , one finally obtains 

^rfo] = S$ c ] + \ Trlog (Gq 1 ) + \ Trlog^o + \ Trlog (l + A F U"Q C )) . (47) 
Eq. (47) can be further simplified to 

13 r$ ] = S$ c ] + \ Tr log (G^\4>S) , (48) 
where, using that Tr log (G^ 1 ) + Trlog(? = Tr log A^ 1 (see Appendix B), 

G^]" 1 = A^ + W'ttc) . (49) 

Notice that the dependence of T on the physical field </>o comes from the nontrivial functional dependence of the 
classical configuration (f> c on <fi . 

V. RENORMALIZATION PROCEDURE FOR THE SINGLE- WELL QUARTIC INTERACTION (D=4) 

The effective action in Eq. (47) is written in terms of non-renormalized parameters. We conclude from the above 
calculation that T has the usual zero-temperature expression when viewed as a function of the classical field <j> c . This 
fact suggests that the renormalization procedure should follow the standard recipe. In this section, we verify that for 
the single-well quartic interaction in d = 3 spatial dimensions. 

We start by introducing a cutoff A for integrations over momenta and adding to T couterterms in order to obtain 
a renormalized effective action, 

(3T R [J }^(3T[$ }-^ J^ x ) E f c{x )-9l J {d ^ x)E ^ c{x) ^9l J { ^ x)e (d^ c ) 2 (x). (50) 

Suitable choices of C\, G2 and G3 should render Tr finite. Renormalization conditions can be fixed using correlation 
functions of the <f) c fields. Each choice of counterterms under that recipe will correspond to a definite set of renor- 
malization conditions in terms of physical correlations involving <p . That translation demands, in principle, the full 
knowledge of <j) c as a function of <p . 
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Let us call vo(x) the saddle-point boundary configuration corresponding to j(x) = 0. Functional derivatives of 
with respect to <j> (x) evaluated at vq(x) lead to the n-point IP I renormalized vertex function. For single- well 
potentials, one finds vq(x) = 0, and r^f^g] admits the expansion 



- °° 1 [ 

t r[4>o\ = / r^ } (a;i,...,a;„)(/) (a;i) ... <f> (x n ) <Pxi . . . d 3 x n , 

n=l U - J 



where 



r^j (xi j • • ■ j x n ) 



O<p i . . .0<p on 

One can also expand T R [(j) a ] in powers of (t> c [4>o\ (using that <p c [0] = 0), 



°° 1 f ~ 

r fl[^o] = / r fl ) ( a; i'---' a; n)^c( a; i) ■•• ^ c (a;™)rf 4 a:i • • • d 4 £„ 

n=l 



where 



5<j) c i ...54> c 



A standard calculation yields 



and 



2 



where 

(2tt) 4 

For massive theories, we use the following zero-temperature renormalization conditions, 



h = J j^r I [A F {k)A F (k + p 3 +P4) + A F (k)A F (k + p 1 +p i ) + A F (k)A F (k + P 3+Pi)} 



1 



- — " w - r ( B } (o,o) = to 2 



(2) 



(2tt) 4 1 5( 4 )(0) dp? 



■(0,0) = 1 



(2tt) 4 1 5( 4 )(0) 



rg ) (o,...,o) = A fi . 



In this case, the counterterms are 

r-A J4 



and C3 = , 



where 



A° F (fc) 



/c 2 + m,Q 



(51) 



(52) 



(53) 

(54) 

(55) 
(56) 

(57) 

(58) 

(59) 
(60) 

(61) 
(62) 



is the zero-temperature free propagator in 4-dimensional Euclidean Fourier space. 

Alternatively, the renormalization conditions can be expressed in terms of correlations of the fields <fi . In principle, 

the general n-point function can be obtained from (53) if one writes 4> c [<i>o] in a power series of (f> , as given by 
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the implicit relation (21). In particular, is a complicated combination of f ^ and On the other hand, the 

2-point function admits a simple expression, 



r^( Pl ,p 2 ) = 0o(Pi) + Q TrA F - d - CsP?) | [^G (r, Pl )]g [3 r ,G (r, - Pl )]g 



(27r)35(p 1+ p 2 ) 

In particular, using the counterterms given in Eq. (61), we obtain, in Fourier space, 



(63) 







(i\ , , Ar tanh /3a;,,, 12 ( /3o;„, 

■ r%\ Pl ,p 2 ) = 2u Pl tanh/?o; Pl /2 + -| — 1 + 



(27r) 3 <5( Pl + ft) * Pl 24/3 V sinh#j Pl . 

When the zero-temperature theory is massless, we adopt the following set of renormalization conditions, 



(2tt) 4 (5( 4 )(0) 



r^( M ,-/i)=^, (64) 



(2tt) 4 ^( 4 )(0) dp\ 

'<*>) 



R 
df {2) 

-^4^,-^ = 1, (65) 



(2tt) 4 (5( 4 )(0) 



r^(Pi,...,P0 



= A K , (66) 



where the arguments of the four-point function are such that J^iPi = 0> Pi = an d (P* + P.?') 2 = 4/i 2 /3 for i ^ j. 
The quantity /i is the renormalization scale. We obtain the counterterms 

Ci = \ f ^ A°(fe) , ^^i/ A (^ A ^ + ri and C3 = - (67) 
The renormalized coupling constant Ap satisfies the following renormalization group equation, 



As a consequence, \r runs with the renormalization scale jj, as 

16tt 2 " & 2n/(3 



A*(/*) = X R (2n/(3) + iXl ^ ] log -4- + 0(X%) . (69) 



(2) 

The function r R ! in the massless case is given by 

f* A > _ , ...... , ,n , *«W tanhfilpJ/2 , 0\ Pl \ 



t2*fs( Pl+Ps > T * ^*rt = *l-*l^ W^RT- l 1 + »J ■ (70) 

The second term on the r.h.s. of (70) plays the role of a mass term for the boundary theory (the first term is the 
kinetic part of that theory). The corresponding dynamical mass comes from the interaction of the boundary field 
with the 77 fields and strongly depends on the external momentum, as shown in Fig. 1. 

One can think of an effective mass of the boundary field as the zero external momentum limit of the two-point 
function, 



This effective mass coincides with the first-order perturbative result for the thermal mass (topt)- The inverse of 

can be thought as the 3-dimensional 
asymptotic behavior of G^ 2 '(|a;i — x 2 \) as 



T^ 2 ' can be thought as the 3-dimensional propagator G^ 2 '(p 1 ,p 2 ) of the field </>q. Numerical calculation shows the 



G ( 2 ) (W)0C ?^P , (72) 



which reinforces the role played by mpr as the mass scale for the boundary theory. 
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FIG. 1: Plot of the dynamical mass (m) of the boundary field <f>o (normalized by the first-order perturbative thermal mass 
mpT = a/Ah/(24/3 2 )) as a function of \p x \ /3, where p 1 is the external momentum. Notice that m approaches mpT in the zero 
momentum limit. 



Naive perturbation theory applied to finite-temperature field theory is doomed to failure due to severe bosonic 
infrared divergences [1]. These are related to zero modes in the Matsubara frequencies, and introduce new mass 
scales associated with collective modes of the system. Incorporating this fact in the description of a thermal system 
in equilibrium naturally leads to the possibility of reorganizing the perturbative series, to different optimizations in 
the diagrammatic formulation, and ultimately to the construction of effective field theories [17, 18]. The process will 
be successful in the measure that one finds appropriate quasiparticles or degrees of freedom to anchor the effective 
theory. 

In this paper we have built an effective field theory that describes the quantum statistical correlations in thermal 
field theories using the expectation value of the boundary fields, which are the effective degrees of freedom, our natural 
stochastic variables and zero modes. This boundary effective theory (BET) is dimensionally-reduced by construction, 
as well as the thermal correlation functions it provides, which can be directly associated with the equilibrium thermal 
system. 

To illustrate our formulation of the BET, we considered a scalar self-interacting field to one loop, but the procedure 
can in principle be extended to gauge fields and we believe that the result holds for higher loop orders. In the scalar 
case, we proved that the one-loop effective action at finite temperature has the same renormalized expression as at 
zero temperature provided it is written in terms of 4> c and zero-temperature propagators are substituted by their 
thermal versions. The classical solution (j) c (r,x) is a functional of the boundary configuration <fio(x) and carries, 
in the underlying statistical mechanical problem, the dynamical information associated to the quantum functional 
integral formalism. Using that connection, one can compute directly all the renormalized thermal vertex functions 
and read off the renormalization group running of all parameters. 

Even at one-loop order, the functional dependence of cj> c on 4>q led to a highly non-perturbativc effective action for 
the boundary field. Our result for T R shows that, although the dynamical mass of the new degrees of freedom the 
boundary fields — have a complicated momentum dependence, its zero external momentum limit coincides with the 
first-order perturbative result for the thermal mass and that has the correct asymptotic behavior as a function 
of |x|, an effect of dynamical screening. These are non-obvious consistency checks of the procedure. 

The rich structure encoded in the d-dimensional vertex functions obtained from the new degrees of freedom of BET 
will, of course, affect the thermodynamic functions and seem to provide a more direct and natural way to compute 
quantities such as the pressure. These results will be presented in a future publication [19]. 
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Appendix A 

In this Appendix, we calculate the derivative of the functional <p c [4>o\- First of all, we multiply Eq. (22) by the 
classical field (P c (t',x'), and integrate over r' and x' . This gives 



c (t, x) = £ dr' J d d x' <j> c { T ', x 1 ) p E + ml] G (r, x; r', x') . 



(Al) 



Now, multiply the equation of motion (Eq. (13)) for (/) c (t',x') by the Green's function Gq(t, x; t', x') and integrate 
over t' and x': 



° = 1 ^ I dVG °( T ' a; ; T '' a; ') [(n E +m 2 )<j> c (T,x') + U'(<j> c (T',x'))] 
Subtracting (A2) from (Al) leads to 



(A2) 



z(T,x) 



dr' / d d x' G (t,x;t',x') 



a E - a E 



~Xt',x') 



- dr' J d d x' G (r, x; r', x') U' 



(A3) 



where the arrows on the differential operators on the first line indicate on which side they act. The first line can be 
rewritten as a boundary term, by noting that 



d 2 - d 2 



B = d»\A 



D 



(A4) 



In eq. (A3), the boundary in the spatial directions does not contribute to the classical field at the point x because the 
free propagator decreases fast enough when the spatial separation increases. Thus, we are left with only a contribution 
from the boundaries in time. At this point, since the boundary conditions for <j) c consist in specifying the value of the 
field at t' = 0, /3, while its first time derivative is not constrained, it is very natural to choose a Green's function G 
that obeys the following conditions, 



G (t,x;0,x') = G (t,x;P,x') = . 
With this choice of the propagator, we obtain formula (21) for c (t, x). From (21), we calculate 
<50 c (t, x) 



(A5) 



J d d x'S(x-y) [d T ,G (x-y,x')i 



- {d u x)' E G {x,x')U"{Ux')) 



,^5cj> c {T',x") 



Ho(y) 



[d T ,G (x;r',y)f - J {d D x)' E G Q {x,x')U" {4>c{x')) 



,(r',x") 

'My) 



(A6) 
(A7) 



Multiplying the previous equation by G 1 (z,x) and integrating over x, we obtain 
(d D x) E G- 1 (z,x) 5 ^2 y ^ = j{d D x) E G- 1 (z 1 x) [d T ,G (x;T',y)f 



- J (d D x) E {d D x)' E G^(z,x)G (x,x')U"(Mx')) Hc }l'; X " ] ■ (AS) 
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Using that Gq 1 (z,x) = 5^{z — x) (D E + m^), we can write 

{UE+ml) ^) +U " {Uz)) ^§) = /(^^M [d T >G ( X ,T>;y)f . (A9) 
From (19), we have 

G- 1 [^ c ](z',z)-^(z'-z)G z - 1 [0 c ] , (A10) 
where G" 1 ^] = D E + m 2 + U"(<p c (z)). Therefore, 

G-'m ^| = J (d D x) E Go\z,x) [d T ,G (x,T';y)f . (All) 

Denoting y by (i"',y) and z by (r",z), notice that 

T< J{d D x) B G^{z,x)G {x;r',y) = d Tl [6^(z-y)6(r" -r')]^ , (A12) 



so that 

G^fec] = *< 3 >(z - y) [d T ,8{r" - r't . (A13) 

Multiplying both sides by G[4> c ](x, z') 5^ (z' — z) and integrating over (d D z) E (d D z)' , we obtain 

^-/(^.(^^(..^W^.)^ (A14) 
= /" dT"d d zG[(b c ](x;T",z)5^(z-y) [d T ,8(r" - r')f (A15) 
= J dr" G[4> c ](x; t", y) [d T ,6(r" - r')f . (A16) 
Integrating by parts, Eq. (23) follows. 

Appendix B 

The free propagator (solution of Eq. (22)) can be explicitly calculated. It is given in Fourier space by 

r ,,,, sinh[cj fc (r> -(3)} sinh(o; fc r < ) 

G (r,T,fc) = — — — , (Bl) 

oj k sinh(/?w fe ) 

where t>(t<) = max (min) {t, t'}. From Eq. (Bl) we obtain 

If, r (rr h y\0 cosh[^ fc (/3/2 - t)] 

[d T ,G (r,r ,fc)] - cosh(wfc/3/2) (B2) 

and 

Q (t, r'; k) = [d T d T ,G (T, t'; fc)]o = 2w fc tanh(w fc 0/2) . (B3) 
It is a simple matter to verify that 



2cj 



sinh[cj fc (T> - (3)] sinh(cj fe T < ) 
uj k sinh(/3cj fc ) 



coshK(/3/2 - t)] coth^ fc /3/2 cosh[u k (/3/2 - t')] 

cosh(cj fe /3/2) 2w fc cosh(w fc /3/2) ' 1 j 



which is the expression in three dimensional Fourier space for 

A F = G Q + d T Go Go l d T ,G . 

It is possible to show that 



ry, , n -\ f d d k sinh/3o; fe 
TrlogG/ = / — 3 log- 



(27T) 2 



Using (B3) and the trigonometric identity 



2tanh- sinhx = e x (1 - e a: ) 2 , 
2 v / j 



we obtain 



^logG^TrlogC^ 



+ log(l-e-^ fc ) 



- TrlogA^+G, 
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(B5) 
(B6) 

(B7) 

(B8) 
(B9) 



where C is an unimportant infinity constant. 
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